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■ Abstract 

The Reggeon field theory with zero transverse dimensions is studied in the Hamilto- 
nian formulation for both sub-and supercritical pomeron. Mathematical aspects of the 
\ model, in particular the scalar products in the space of quantum states, are discussed. 

Relation to reaction-diffusion processes is derived in absence of pomeron merging. Nu- 
merical calculations for different parameters of the models, a(0) — \ = [i and the triple 
pomeron coupling constant A, show that the triple pomeron interaction always makes 
amplitudes fall with rapidity irrespective of the value of the intercept. The smaller the 
values of the ratio \j \x the higher are rapidities y at which this fall starts, so that at 
I small values of A it begins at asymptotically high rapidities (for A//i < 1/4 the fall is 

noticeable only at [iy > 100). No visible singularity is seen for the critical pomeron. A 
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perturbative treatment is proposed which may be useful for more realistic models. 

1 Introduction 



X 

At high energies the Quantum Chromodynamics tells us that particle interaction may be 
mediated by the exchange of hard pomerons, which are non-local entities propagating ac- 
cording to the BFKL equation and splitting into two or merging from two to one with the 
known triple pomeron verteces. This most complicated picture can be drastically simplified 
in the quasi-classical treatment supposedly valid when at least one of the colliding hadrons is 
a heavy nucleus. This leaves only tree diagrams which can be summed by comparatively sim- 
ple equations, which at least admit numerical solution. However at present much attention 
is drawn to the contribution of loops. This is obviously a problem of a completely different 
level of complexity, since it amounts to solving a full-fledged non-local quantum field theory. 
At present a solution of even some local quantum field theory has not been known except 
for the one-dimensional case. Naturally for the study of loops in QCD an one- dimensional 
model has been lately chosen as a starting point. 

The one-dimensional ("toy") model for the QCD interaction was proposed and solved 
by A.H.Mueller [lj. However much earlier a model similar in structure was studied in 
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the framework of a local Reggeon field theory (RFT) with a supercritical pomeron and 
imaginary triple pomeron interaction in a whole series of papers [21 El SI EJ, [T| As expected 
the one-dimensional reduction of the RFT admitted more or less explicit solution and lead 
to definite predictions as to the behaviour of amplitudes both for hA and AB interaction 
(of course with the understanding of a highly artificial physical picture, very far from the 
realistic world). In the course of study many ingeneous tricks and techniques were employed 
and the underlying complicated dynamics was revealed. However some subtle points have 
remained unclear in our opinion. In particular the nature of correspondence between the 
functional and Hamiltonial approaches and the origin of restrictions on the form of wave 
function and the range of its variable have not been fully exposed. 

To elucidate these points is the primary aim of this note. We also try to follow the line 
presented in [HI [9] to relate the RFT with the probabilistic parton picture and the so-called 
react ion- diffusion model for the case when there is no merging of pomerons (fan diagrams). 
This gives us a possibility to find the contribution from the loops with the highest extention 
in rapidity just by joining two fans propagating from the projectile and target towards the 
center. We additionally develop a perturbation theory for small A to compute the evolution 
by noting the PT-symmetry of the model. We finally report on some numerical results 
on solving the RFT with all loops, which can be obtained quite easily by integrating the 
corresponding evolution equations. 

Note that we study exclusively a model with the triple pomeron interaction only. A 
different model which additionally includes a quartic pomeron interaction with a particular 
fine-tuned value of the coupling constant has lately received enough attention in literature 
[HI EE HOI [11] . In our opinion the quartic interaction has no immediate counterpart in the 
realistic pomeron theory with a large number of colours N c , where it is damped by factor 



2 Functional integral formulation 

The toy model of Pomeron interaction with a zero- dimension transverse space can be defined 
by a generating functional 

Z = J D(f)D^e- s } (1) 

where 

S = J dy£, (2) 

and 

£ = i(0V y - 0t0) _ + *A0V(0t + 0) _ i(Jo t + jy 0) . (3) 

Here [i is the pomeron intercept (a(0) — 1). For the supercritical pomeron > 0. Triple 
pomeron coupling constant A is also positive. The sources j = $(y)9i am d jy = 5{y — Y)g 2 
represent the interaction with the target and projectile respectively. In the real world g\ = 
giip) = AgT A {b) and g 2 = #2(6) = BgT B {b) where g is the pomeron- nucleon coupling 
constant, A and B are atomic numbers of the colliding nuclei and T4 and Tb their profile 
functions. 
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The functional integral ([I]) converges for // < provided the two field variables and 0t 
are complex conjugate to one another. In fact putting = 0i + 202 and 0^ = 0i — 202, with 
both 0i 2 real, we find that all terms in the Lagrangian ([3]) are pure imaginary except for 
the mass term proportional to 

0V = 0? + 02 > 0. 

So the integrand in ([I]) contains a factor exp (ji J dy{(j)\ + 0|))> which guarantees convergence 
if fi < and the integration over the fields 0i and 02 goes over the real axis. For /i > the 
integral does not exist. Note that the integral does not exist for any sign of /i either, if one 
changes the integration contour in 0! and 02 and so over and 0t. This means that it is 
impossible to pass to new fields u = i<$ and v = i<p and take u and v real directly in the 
functional integral as in [3], since it requires unlawful rotation of the integration contour in 
the integration over 0i. 

Obviously the functional formulation can serve only to define the model for the subcritical 
pomeron with /i < 0. To define the theory for positive values of \i one has to recur to analytic 
continuation in fi. Note that the properties of the functional integral make one think that 
there is a singularity at /i = 0. In a fully dimensional RFT it was conjectured that this 
singularity was related to a phase transition. However already in [6] it was argued that there 
was no phase transition in the zero-dimensional world and the amplitudes were analytic in 
/i on the whole real axis. As we shall see by direct numerical calculation, there is indeed no 
singularity at /i = in scattering amplitides. 

In fact the functional formulation is only supported by the fact that it reproduces the 
perturbative diagrams for the pomeron propagation and interaction. In the following section 
we introduce an alternative, Hamiltonian formalism, which gives rise to the same perturba- 
tive diagrams (such an approach as a starting point was briefly mentioned in [1]). However 
in contrast to the functional approach it does not involve any limitation on the value or sign 
of the intercept fi. Since for \i < the perturbative series seems to be convergent, the two 
formulations are completely equivalent for the subcritical pomeron. Therefore the Hamilto- 
nian formulation gives the desired analytic continuation of the model to positive values of 
/i. 



3 Hamiltonian formulation 



3.1 Basic definitions 

Forgetting the functional integral ([1]) altogether, we now start from a quasi-Schroedinger 
equation in rapidity: 



dy 

and postulate the form of the Hamiltonian H 



-H*(y) (4) 



H = -^0 + ^(0 + 0^)0 (5) 
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as a function of two Hermithean conjugate operators <fi and $ with the commutator 

[0,01 = 1- (6) 

The Hamiltonian is obviously non-Hermithean. We standardly split it into a free and inter- 
action parts: 

H = H + Hj, # o = -/i0 f 0, H I = i\^{<f> + ^)<f>. (7) 

Next step is to define physical observables in this picture. In accordance with the com- 
mutation relation ((J6]) we define the vacuum state normalized to unity, by the condition 

0^o = 0. (8) 

All other states will be built from \l/ by application of some number of operators 0^. We 
postulate that the transition amplitude from the initial state at rapidity y = to the 
final state at rapidity y is given by 

iA fi = (*/|*i(y)> = (*/|e-*»|¥ 4 ). (9) 

Here 

ttife) = e-^tti (10) 

is the initial state evolved to rapidity y. The amplitude Afi is imaginary positive so that the 
matrix element on the right-hand side of (9) is negative. 

It will be convenient to define the initial and final states in terms of some operators acting 
on the vacuum state: 

^ = ^(0^0, * f = F f (tf)* . (11) 
This allows to rewrite the amplitude as a vacuum average 

iA fi = (F;(0)e-^(0t)) ; (12) 
where we do not explicitly write the vacuum state m which the average is taken. 

3.2 Perturbative expansion 

We are going to prove that the perturbative expansion of this theory leads to the same 
Feynman diagrams as from the functional approach. To this end we present 

iAfi = (e Hoy F}((j))e- Hoy e Hoy e- Hy F0)). (13) 

From it follows that 

<j>(y) = e Hov <pe- Hoy = e m <p. (14) 

The operator 

U(y) = e Hoy e- Hy (15) 

satisfies the equation 

= e Ho V{Ho _ H)e -m = _ e HoV H . e -H y u{y) = fy^y) (1 6 ) 

dy 
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where Hz(y) is the interaction in the interaction representation: 

Hi{y) = e^Hj^, ^)e~ H ™ = H^y), 0%)), (17) 

with the rapidity dependent operators (14) and 

<j)\y) = e Hoy ^e- Hoy = e'Wtf. (18) 

Note that at y > the operators <fi(y) and (^{y) are no more Hermithean conjugate to 
each other, since the normal time is changed to the imaginary one. However the similarity 
transformations (fl4l) and (fl8l) do not change the commutation relations. The solution of 
ffT6l) with the obvious initial condition U(0) = 1 is standard 

U(y)=T y ex P {-J Q V dy'H l (y')}, (19) 

where T y means ordering the operators from left to right according to decreasing y's. So we 
get the expression for the amplitude as 

iAfi = (F;(e^0)T y exp { - jT dy'H^F^)). (20) 

Expansion of the T y exp gives rise to standard Feynman diagrams. Verteces for pomeron 
splitting and merging will be given by the expected factors iX. The pomeron propagator will 
be given by 

P( Vl - y 2 ) = (Tytyfa), <f>\y 2 )}) = efto^gfa _ y2 ), ( 2 1) 

which is the correct pomeron propagator in our toy model. So indeed the theory defined by 
this Hamiltonial formulation gives rise to standard Feynman diagrams for the RFT. 

Note that the Hamiltonian formulation is based on the evolution equation in rapidity 
and does not require the intercept fi to have a definite sign. It looks equally good both for 
positive and negative values of /i. For negative \i it produces a perturbative expansion which 
is identical to that in the functional approach. So it gives the desired analytic continuation 
of the functional approach to the region of positive fi, where the latter looses sense. 



3.3 Passing to a real Hamiltonian 

It is perfectly possible to continue with the originally defined Schroedinger operators <j) and 
0' which are Hermithean conjugate to each other and satisfy the commutation relation ([6]). 
However it is convenient to pass to new operators u and v in terms of which all ingredients 
of the theory become real. Following the old papers (starting from [2]) we define 

u = i<p\ v = i(f). (22) 

The new fields are anti-Hermithean to each other 

= —v, v' = — u (23) 
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and satisfy the commutation relation 

M = -l. (24) 

However v^q = and the states are generated by action of u. So v is the annihilation and 
v is the creation operator with abnormal commutation relation ( [24]) . 
In terms of new fields the Hamiltonian becomes real 

H = fiuv — Xu(u + v)v. (25) 

It is important that also the operators creating the initial and final states become real. It is 
normally assumed that for the scattering of two nuclei the initial and final wave functions 
should correspond to the eikonal picture. In terms of creation operator u this leads to 

* i(f) = F i(f) (u)* Q , F i(/) (u) = l-e-*W\ (26) 

So the expression for the amplitude is rewritten in terms of real quantities: 

iAfi = (F f (u)y \e- H y\F t (u)V ) = (V \F}(-v)e- H yF t (u)y ) 

= <(l - e-^ v )e~ Hy (l - e~ 9iU )). (27) 
Since H^/q = the term independent of $ and gj vanishes, so that we can also write 

iA fi = -(e-'f v e- H »(l - e-*")> = -{^~ afV Hv, «)*o), (28) 
where Fi(y,u) is the operator which creates the evolved initial state. Since 

%(y) = e' Hv % = e- Bv Fi(u)* = F^y, «)#„, (29) 

we have 

^^- = -H{u,v)F i {y,u) (30) 

with an initial condition 

F,(0, M ) = l-e-^. (31) 
The commutation relation (]24j) allows to represent 

v = ~l (32) 

and then ([28]) implies that to find the amplitude one has to substitute u by in Fi(y,u) 

iA }i = -F i {y,g f ). (33) 

This gives a practical recipe for the calculation of the amplitude. One has to solve Eq. 
(1301 with the Hamiltonian (1251) in which v is represented according to ( [32]) : 



(9 d 

H = ~^u— + \u 2 — - \u— (34) 
ou ou Ou z 

and the initial condition f[3T|) . After the solution is found, one has to substitute u by in 
it. 
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3.4 Target-projectile symmetry 

Taking the complex conjugate of (1281) we find 

- iA* fi = ((l - e^)e- H 'y(l - <?*«)) = iA if (X -> -A,^ (/) - (35) 

Having in mind that the amplitude is pure imaginary, we see that interchanging the target 
and projectile leads to changes of sign in the triple pomeron coupling and the couplings to the 
external particles. On the other hand changing u — > —u in Eq. ( 130]) with the Hamiltonian 
( I34p and the initial condition fl3T|) and also changing signs of A and Qi obviously does not 
change the solution: 

Fi(y, -u, -g h -A) = F^y, u, g h A), (36) 

from which we find 

Fi(y, ~9f, -gh - A) = Fi(y, g f , g h A). (37) 
So the interchange of the target and projectile does not change the amplitude. 



3.5 Scalar products 

As follows from our derivation the exact mathematical realization of the scalar product in 
the space of wave functions \1/ is irrelevant for the resulting formulas. The representation of 
the operator v as a minus derivative in u serves only to express its algebraic action on the 
operator u in accordance with the commutation relation. It does not require to introduce 
a representation for the wave function \1/ in which u is represented as multiplication by a 
complex number and correspondingly the scalar product is defined by means of integration 
over the whole complex plane u with a weight factor e~ uu * aas in [6j. One may instead 
represent operators u and v in terms of the standard Hermithean operators p and q as 

u= -j={q-ip), v = -^={q + ip) (38) 

and define a scalar product as 

(*!|tt 2 ) = J dq^l(q)^ 2 (q) (39) 

with p = —id/dq. The vacuum state \l/o wih be given by the oscillator ground state in this 
representation. 

However this is not the end of the story. New aspects arise in the process of solution of 
the Eq. ( 130]) for Fi. The equation itself is regular everywhere except possibly at u = 0, where 
however it is also regular having behaviour as a constant or u. We are obviously interested 
in the last behaviour since from the equation it follows that Fi(y, u = 0) = 0. So in principle 
our solution turns out to be regular in the whole complex w-plane. Formally this means that 
one may choose any continuous interval in this plane and evolve the initial function given in 
this interval up to the desired values of rapidity. The immediate problem is what interval 
one has to choose. The final expression for amplitude requires that the value of gf has to 
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belong to this interval. In particular for positive gj the interval should include at least a 
part of the positive axis. 

Numerical calculations show that it is indeed possible to find a solution to the evolution 
equation provided the initial interval is restricted to only positive values of u (including the 
desired value u = gf). 

To understand this phenomenon one has to recall some results from the earlier studies. It 
was found in [I] that there exists a transformation of both the function and variable which 
converts the non-Hermithean Hamiltonian (j34|) into a Hermithean one. However it is only 
possible for values of u covering the positive axis rather than the whole real axis where the 
equation should hold. Namely one introduces variable x by 

u = x 2 > (40) 

and makes a similarity transformation of the Hamiltonian 

H = W- l HW, W(x) = v^e^ 4 -^ 2 

to find the transformed Hamiltonian as 

where we standardly denote p = fi/X. The Hamiltonian if is a well defined Hamiltonian on 
the whole real axis with a singularity at x — 0. The interaction potential infinitely grows at 
\x\ — > oo. So restricting to x > and requiring eigenfunctions to vanish at the origin and 
at x — > oo one finds a discrete spectrum, starting from a certain ground state Eq to infinity. 
The transformed initial function 

FAy = 0,x) = -Le-^ 4 ^* 2 ( 1 - e"^ 2 ) (43) 

\/x v ' 

is well behaved both at x = and x = oo and at < x < oo can be expanded in the 
complete set of eigenfunctions F n of H: 

Fi{y = 0,x) = ^(Fi{y = 0,x)\F n )F n (x) 

n=0 

Here a new scalar product has been introduced for the transformed functions: 

(iq|F 2 ) = / dxF*(x)F 2 (x) 
Jo 

Evolution immediately gives 

Fi(y,x) = ^(^(y = 0,x)\F n )e~ Eny F n (x) 

n=0 

and applying the operator W one finds the solution as 

F(y,u) = J2(F l (y = 0,x)\F n )e- E ^W(x)F n (x), u = x 2 . (47) 

n=0 
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(41) 



(42) 



(44) 



(45) 



(46) 



Note that we could just as well start from negative values of u and define instead of (HOI) 



u = -x • , x > (4£ 



Doing the same transformation of the Hamiltonian we then obtain (1421) with the opposite 
overall sign and an opposite sign of p, that is the transformed Hamiltonian is now —H(p — > 
— p). All eigenvalues now are negative and go down to — oo. If we take the expansion of 
the initial function analogous to (|4*4l) and try to evolve it we obtain a sum with infinitely 
growing exponentials in y. Obviously such evolution has no sense. This explains why with 
the positive or negative signs of A and g's one can start evolution only from correspondingly 
positive or negative values of u, but never from intervals including both. 

Also note that considered as a function of u, eigenfunctions are analytic in the whole 
w-plane. Starting from, say, u > their values at u < can be obtained by analytic continu- 
ation to pure imaguinary values of variable x. Obviously eigenfunctions should exponentially 
grow as u —>■ — oo (otherwise we shall find positive eigenvalues of the continued Hamiltonian 
—H(p — > —p)). Similarly starting from u < and passing to positive u one finds that eigen- 
functions exponentially grow at large positive u. So considering the initial Hamiltonian H 
on the whole real u axis we find that its spectrum goes from — oo to +oo with eigenfunctions 
which vanish only either at large positive u for positive eigenvalues or at large negative u 
for negative eigenvalues and exponentially grow otherwise. This prevents introduction of the 
Bargmann scalar product with integration over the whole complex w-plane. 

Finally let us investigate the qualitative behaviour of the system in the large A/p = 1/p 
limit, which may be of interest in the generalization to more dimensions, where the effective 
A is (infinitely) large (see [5|). For this purpose the representation given in Eq. (142]) is 
convenient, since it admits perturbative expansion in p. In the limit p = one has 

which means that the amplitude will behave as e~ yE °, where E is the ground state energy 
oiH. 

One can perform a semiclassical estimation of the ground state of the Hamiltonian in Eq. 
( 149]) by imposing the condition 

^i^l^H-),, (50) 

where V(x) = + x 6 — 2x 2 and X{ are the positive real solutions of the equation V(x) = 
4E/X. For n = one obtains E ~ 3.7 A/4 > 0, which is in a very good agreement with the 
results of numerical evolution for X/ p, > 3, which indeed shows a decay of the amplitude as 
e~ yE ° after a short initial period of evolution. 

Note that during this initial period one observes a rise of the amplitude, which may seem 
strange in view of the fact that all eigenvalues of the transformed Hamiltonian ( T4T]) are in 
fact positive. This rise is explained by the interference of different contributions which may 
enter with different signs due to the non-unitarity of the similarity transformation (|4~T]) from 
the initial non-Hermithean Hamiltonian. 
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4 No pomeron fusion: fan diagrams 
4.1 Solution 

A drastic simplification occurs if in the Hamiltonian one drops the term corresponding to 
fusion of pomerons ('fan case'): 

d 

H — > Hf an = fiuv — Xu 2 v = — {pu — Xu 2 ) — . (51) 

Then, as has been known since [12] , the solution can be immediately obtained in an explicit 
form. Indeed, following [8], define a new variable x by 



du 
pu — Xu 2 



d* = .... ,., a . (52) 



which gives (up to an irrelevant constant) 



1 1 U ^ /KQ\ 

z = —ln , p — — . (53) 

[A u — p X 

The inverse relation is 

u = — . (54) 

1 - e-i» 

Our evolution equation for F takes the form 

,d_ d_ 

^dy dz 

with a general solution 



F(y,z) = (55) 



F(y,z) = f(y + z) (56) 
where / is an arbitrary function. It is to be found from the initial condition: 

F(0, z) = f(z) = 1 - e~ 9iU = 1 - e ~ w i=^. (57) 

From this we obtain [S] 



P 



F fan (y, u) = 1 - exp ( - — = 1 - exp [ - i + • (5 



9iue 
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This solution is evidently non-symmetric in the projectile and target. In practice one consid- 
ers such an approximation in the study of hA scattering. Then it represents the sum of fan 
diagrams propagating from the projectile hadron towards the target nucleus. This situation 
corresponds to the lowest order in powers of gi [12] : 



Ffan = - If , (59) 
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We recall that the amplitude is obtained from ([58]) or (1591) by putting u = gf 

We note that the solution (1591) is not analytic in the whole complex u plane: it evidently 

has a pole at some negative value of u, which tends to zero from the negative side as y — ► oo. 

So at least for this solution one cannot formulate a scalar product of the Bargmann type 

with integration over the whole complex u plane. As we have stressed, the existence of such 

a scalar product is not at all needed. 

In fact the Hamiltonian Hf an has the same spectrum as the free Hamiltonian H (Eq. 

([7])). The two Hamiltonians are related by a similarity transformation 

H fan = VH V-\ V = e- u2v 'p. (60) 

It follows that eigenvalues of Hf an are all non-positive: E n = —[in with n = 0,1,2,.... 
However the amplitudes (|B^|) and (IB^I) do not grow but rather tend to a constant at large 
y. As seen from the structure of expressions (l58l) and (l59l) this is a result of interference of 
contributions with different E n , which separately grow at high y. One can reproduce the 
result of the fan evolution using the form in Eq. (1601 by noting that e~ yHa is an operator 
which shifts the variable lnw by /i while V corresponds to a shift in the variable 1/u by 1/p. 

The fan evolution is to be compared with the situation with the total Hamiltonian H, 
whose eigenvalues are all non-negative for the branch with positive u. In this case the 
solutions generally fall at large y and became nearly constant only at very small A when 
one of the excited levels goes down practically to zero due to a particular structure of the 
transformed Hamiltonian (H2l) |4J . So the mechanism for the saturation of the amplitudes 
at large rapidities is completely different for the complete Hamiltonian H and its fan part 

Hfan- 



4.2 Relation to the reaction-diffusion approach 

The fan case admits a reinterpretation in terms of the so-called reaction-diffusion processes 
[9], which may turn out to be helpful for applications. Let us present a solution F(y,u) as 
a power series in u 

F(y,u) = -£c n (y)u n . (61) 

n=l 

The evolution equation gives a system of equations for coefficients c n (y): 

dc ™ V = nncn{y) - \{n - l)cn-i(y)- (62) 
Consider first the case \i > 0. We rescale y = yj \i and present 

c n {v) = -A~pf- n v n {y) (63) 

to obtain an equation for v n 

= nVn{ y) + n{n _ l)u n ^(y). (64) 
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This equation is identical with the one which is obtained for multiple moments of the prob- 
ability Pk{y) to have exactly k pomerons 



^(y) = £n(y)7-r— v ( 65 ) 

k=n v i' 



provided the probabilities obey the equation 

dP n (y) 



dy 



-nP n (y) + (n-l)P n ^(y). (66) 



The equations for P n are such that they conserve the total probabilty assumed to be equal 
to unity. Also if the initial probabilities are non-negative, they will remain such during the 
evolution. This gives the justification for the interpretation of P n as probabilities. 

For further purposes it is convenient to introduce a generating functional for the proba- 
bilities P n 

Z(y,u)=Y,Pn(y)u n . (67) 



In terms of this functional 



1 / d \ n / d \ n , 

p « = ^Gr z > v »={jr ) z ■ 68 

n\ K ou J u=o K ou J u=i 

From the system of equation for P n one easily finds an equation for Z(y,u) 

dZ(y,u) 2 ^ dZ(y,u) 

^y— = {U - U) ^^- (69) 

Note that this is essentially the same equation as for F with parameters put to unity and 
the opposite sign of the Hamiltonian. So its solution is immediately obtained from fl56|) by 
changing sign of y 

Z(y,u) = z(0,z(u)-y), u = ^—^ (70) 
For the simple case of hA scattering this gives 



7hA 



(71) 



l + u{e~y - 1) 
Expansion in powers of u gives the probabilities P n as 

p% A = e -y(i - e -ff) n_1 . (72) 

For a general case we first construct Z(y = 0) as a sum 

in — l) n 

Z(0,«) = l + 5>»(°)- ^- (73) 

n=l n - 
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with values of z/ n (0) known from the initial distribution: 

u n (0) = p^gf. (74) 
Subsequent summation and shift z — > z — y lead to the final result 

11 / (u— l)e s \ 

Z(y,u) = 1 b-explgip — -). (75) 

p p v u — (u — l)e yJ 

The probabilities should be obtained by developing this expression around u = 0, which is 
not easily done in the general form. 

Now let us briefly comment on the case of the subcritical pomeron p = — e < 0. We again 
rescale 

y = y'/e (76) 

and instead of f[6"3"j) put 

cn(y r ) = -/- n K(y') (77) 

n! 

where now p = — e/A < 0. The resulting equations for u' n are 

dv' (v') 

= -nK(y') + n(n - iV^y'). (78) 

As we see they become similar to the equations for P n in the case of p > and in fact coincide 
with them for the quantities v n = v' n jn\. This means that for negative p the relation between 
the RFT and reaction-diffusion processes is different from the case of positive p. In fact the 
interpretation of the amplitude in terms of P n and v n interchanges: the rescaled coefficients 
c n in the amplitude give directly the probabilities to find a given number of pomerons and 
their multiple moments have to be defined in terms of these quantities in the way similar to 
the construction of probabilities from their moments for p > 0. 



5 Big loops 

One can use the Hamiltonian approach to find the contribution of loops with the maximal 
extention in rapidity, similar to the approach first taken in [13] in the fully dimensional 
dipole picture. This issue has been also investigated for a different model with a specific 
quartic interaction where an exact solution can be built [10J. To this end we first split the 
evolution in two parts: 

iA fi {y) = {F}(-v)e- H yF^u)) 

= ( J P;(- W )e-^/ 2 e-^/ 2 F i ( W )) = ((e-^ 2 F f (u)^e- H ^ 2 Fi(u)). (79) 

Next we neglect fusion of pomerons during the first part of the evolution and merging of 
them during the second: 

iA fi (y) = ((e- H )- 2//2 J P / («)) t e-^» 2//2 F i (w)). (80) 
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Here Hf an is given by Eq. (l5TT) and H\ an differs from it by the sign of A. This approximation 
takes into account loops which are obtained by joining at center rapidity the two sets of fan 
diagrams going from the target and projectile. 

Both operators inside the vacuum average are known and given by Eqs. ( 1581) or (1591) 



with the understanding that for the solution with Hj an one has to change the sign of A. To 
simplify we restrict to the lowest order in both coupling constants (?«(/)• Separating g^gj we 
then obtain the pomeron Green function with loops of the maximal extention. In this case 
the right factor in flHUl) is given by 

e -H fan y/2 u = «f 

i + Hr e w/2_i) 

and the right factor is 

(e f anyi u) = 7 r- = 7 r-. (82) 

So the pomeron Green function is 

G(y) = -( 4- V -r \)- (83) 

Its calculation is trivial. Denoting 

6= i(>//2_f) ( 84 ) 

we present the left operator as 

VCml2 e^^— = e^ 1 -(l- 1 -^-). (85) 



1 + p( eA 



MS//2 _ 1 ] 1 + bv 6 V + l 



Action of the unit operator gives zero so that we are left with 
Now we present 

I nno 

dxe -x(v+l/b) £ 87 J 



v+l/b Jo 

and, since operator e~ xv just substitutes u by x in the vacuum average, obtain 

1 r°° , _,/!. x ,„, 1 /. 1 



G (,) = e «- f rfxe^^- = + ^El (-1, 



!-,2- 



In the asymptotic limit y — > oo, in this approximation, the propagator tends to a constant 
value 

G(y)y-+oo ^ P 2 — T2 ^ P 2 - (89) 



(e«/ 2 - 1 
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The loops tamed its initial exponential growth as e m illustrating the saturation phenomenon 
in the old Gribov's RFT with the supercritical pomeron. 

Note that the same result can be obtained by developing in the number of interactions [HI 
Ej. We find worth to remind to the reader how one proceeds in this case. In fact we can 
expand both wave funcions in powers of u and v respectively as in (I6T1) . Obviously only 
terms with equal numbers of u and v contribute with 

(v n u n ) = (-l) n n\. (90) 

So we get 

G(y) = -J2(-l) n nlcl(y/2), (91) 

n=l 

where from (1591) we find 

c n (y/2) = e ^ 2 {-b) n - x . (92) 

We find a divergent series 

G(y) = -^y^(-b 2 ) n n\ (93) 

n=l 

which is however Borel summable. Presenting 

POO 

n\ = / dte~H n (94) 
J o 

we get 

I poo POO f 

G(v) = - r «l *^ £(-*)• = e-j[ dte-'^. (95) 

This is the same expression (1881) which we obtained before. So these simple but not very 
rigorous manipulations give the correct answer. 

It is instructive to express the same result in terms of the probabilities P n . Actually it 
is a straightforward task. Rescaling the coefficients c n in (1911) we rewrite this expression in 
terms of multiple probability moments v n 

^)=vi:(4)^/4' ( % ) 



n=l 



Now we express the multiple moments in terms of the probabilities according to (|66l) to 
finally obtain 

G(y) = J2 PkivmPtvmFu, (97) 
k,i=i 

where F^i are coefficients independent of energy 

F„ = X)n!Cycr(-i) B " 1 . (98) 

n=l P 

It is easy to find the asymptotical values of Fki at large k and I. Then we can approximate 

Ck-^, n « k -> oo. (99) 
15 



Under this approximation 

^v^(-S"=/(i-.-'*) (ioo) 

n=l y 

Then we get from fl96|) 

G(y) = P 2 (l - E P fe (y/2)P,(l//2)e^ 2 ) (101) 
fc,/=i 

At high y and fixed n the probabilities P n (y) become independent of n and equal to e _/4/ 
(see Eq. ( !72|) ). Then we get 

G{y) = p 2 (l-e-M E e- fc ^ 2 ), (102) 



k,l=l 



with the correct limiting value p 2 at ?/ — > oo. 



6 Perturbative analysis for small A and PT-symmetry 

Having in mind possible generalizations to the realistic two-dimensional world, in this section 
we study perturbative treatment of the model at small values of A. As follows from the 
previous studies [U, El [6] actually the point A = corresponds to an essential singularity in 
the Hamiltonian spectrum corresponding to the existence of a low- lying excited state with 
a positive energy 

AE = J^e-fl 2 . (103) 

This states dominates evolution at very high rapidities. However at earlier stages of evolution 
this unperturbative component may play a secondary role as compared to purely perturbative 
contributions. We also point out in this section some general properties of the Hamiltonian 
H given by Eq. (5), which belongs to a class of non-Hermithean Hamiltonians widely studied 
in literature. 

Inspecting the Hamiltonian H we can see that even if it is not Hermithean, it describes 
a PT-symmetric system evolving in the imaginary time. This property has been intensively 
studied mainly for anharmonic oscillators |15j . Indeed H is invariant under the product of 
parity and "time" -reversal transformations: — > — 0, (ft — > and i — > —i. 

As the transformation to the form (42) has shown, the spectrum of H is real. This 
fact is true whenever for a diagonalizable operator H certain conditions are fullfilled [16J. 
The most important is that H should be Hermithean with respect to some positive-definite 
scalar product (-|-)^ = (-|^-), constructed using a positive-defined metric operator r\ = e~ Q 
with Hermithean Q. This defines the ?7-pseudo-Hermiticity property of H, which may be 
conveniently written as 

H ] = r] H 7]- 1 = e- Q H e Q . (104) 

Considering the set all positive-definite metric operators, it is easy to see that the parity 
operator P belongs to it. Then one can define a C operator by C = P _1 r/ = r/ _1 P which 
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commutes with the Hamiltonian H and the PT symmetry generator. In this formalism one 
should therefore introduce as a scalar product the CPT-inner scalar product 



(105) 



The r] operator is in general unique up to symmetries of H. It describes the nature of 
the physical Hilbert space, the observables being ?7-pseudo-Hermithean operators. Using a 
symmetric (second) form of the scalar product (105) we define 



h = e- Q ' 2 He Q ' 2 



(106) 



which is Hermithean with respect to the scalar product in this form. The eigenvalues of h 
and H obviously coincide. 

The operator Q can be determined perturbatively. In full generality we write 



H = H + A#i 



(107) 



and in the the perturbative expansion Q = XQi + \ 3 Qs + • ■ ■ we may derive the terms Qi 
using Eq. (I104p and matching the terms at each order in A. Here we shall restrict to the 
first two orders in perturbation theory 



°' ~2~ 
H , Y 



1 

~3 



2 . 



Qi 

2 



(108) 



Inserting this into (106) for for h = + X 2 h^ + \ A h^ + ■ • ■ one then obtains 





= #o 






1 

~ 2 


Hi, 


Qi 

2 




1 

~ 2 


Hi, 


Qz 

2 



TT Qz 

fin, — 



Qi 

2 



(109) 



Let us now use these results for the system under investigation. We define the operators 



and write Hn 



N = (j) ] (t) , R n = N n 4> + ^N n , S n = N n (f)-^N r 
—fiN and Hi = iR\. One easily finds that 



Qi 

2 



Si 



Qz 
2 



S 2 + S! 

3 



Inserting this expressions into Eq. (11091) one finally obtains for h up to order A 4 : 



(110) 



111) 



h 



-H \ N + * (3N 2 -N) + -^ (-12N 3 + QN 2 - 2N + Ah^) + O (^) 



(112) 
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where 




(113) 



is a term not diagonal in the number operator N and thus contributes to the energy eigen- 
values only at a higher order in A///. 

Note that the transformation (|41|) for positive u and A leads to Hamiltonian (J52J), which 
has all eigenvalues positive. In contrast the perturbative eigenvalues are seemingly non- 
positive for small enough A and certainly such at A = 0. This shows that the perturbative 
series is not convergent. Note however that it provides a reasonable asymptotic expansion 
at small values of A for not very high rapidities and describes well the initial period of evo- 
lution when the amplitudes rise. As an illustration, using this perturbatively approximated 
spectrum one may compute the pomeron propagator at two loops using 



expanded in the basis \n)h of the eigenstates of h {h\n) h = E n \n) h ). Comparison to the 
exact numerical evolution in the next section shows a perfect agreement up to some value 
of fiy, much below the saturation region, beyond which the perturbative approximation, 
which continues to grow in rapidity, fails. In particular for X/fi = 1/10 the perturbative 
approximation is very good up to fxy < 3. 

In the limit A — > one finds that the initial period of evolution, when the amplitude 
grows, extends to infinitely high rapidities and the propagator continuosly passes into the 
free one, growing as exp(fiy). So in spite of the essential singularity at A = 0, the amplitudes 
seem to have a well defined limit as A — > 0, when they continuously pass into free ones. 

One might expect reasonable to apply a similar perturbative approach for some interval 
of rapidity also in the case for non zero transverse dimensions. 

7 Numerical calculations 

As mentioned in the introduction, most of the qualitative results in the toy model were 
obtained thirty years ago [21 El lU El El E] • They were restricted to the case of very small 
triple pomeron coupling A, when the dynamics becomes especially transparent and admits 
asymptotic estimates. However the smallness of A in the zero-dimensional case actually does 
not correspond to the physical situation in the world with two transverse dimensions. If 
the transversce space is approximated by a two-dimensional lattice with intersite distance a 
then the effective coupling constant in the zero- dimensional world is inversely proportional 
to a and thus goes to infinity in the continuum space (see [5]). So the behaviour of the 
model at any values of \x and A has a certain interest. In particular it may be interesting 
to know the behaviour at \i — > 0, since the functional integral defining the theory for fi < 
diverges at // > and one can expect a certain singularity at \i = 0. The validity of different 
approximate methods to find the solution may also be of interest in view of applications to 
more realistic models. 

Modern calculational facilities allow to find solutions for the model at any values of \i 
and A without any difficulty Our approach was to evolve the wave function in rapidity 



G(y) = (O^e-^IO) = (0\<f> e Q/ V^e- Q/2 f |0> , 



(114) 
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directly from its initial value at y = using Eq. (I30I) by the Runge-Cutta technique. This 
straightforward approach proved to be simple and powerful enough to give reliable results 
for a very short time. True, the step in rapidity had to be chosen rather small and correlated 
with the step in u. In practice we chose the initial interval of u as 

0<m<20 (115) 

divided in 2000 points. The corresponding step in rapidity had to be taken not greater than 
2.5 ■ 10~ 3 

Presenting our results we have in mind that amplitudes depend not only on the dynamics 
but also on the two coupling constants g^fi . To economize we therefore restrict ourselves to 
show either the full pomeron propagator 

P( y ) = -(ve- Hv u) (116) 

or the amplitude for the process which mimics hadron-nucleus scattering 

A(y,g) = -(e-^e- H y U ) (117) 

(with gj = g > and g>j = 1). 

7.1 Pomeron propagator 

Straightforward evolution according to Eq. (|30|) from the initial wave function Fi(0,u) = u 
gives the results for the pomeron propagator shown in Fig. 1 for different sets of p, and A. 
We recall that the free propagator is just 

Po(y) = e™. (118) 

Inspection of curves in Fig. 1 allows to make the following conclusions. 

1) Inclusion of the triple-pomeron interaction of any strength makes the propagator fall 
with rapidity. This fall is the stronger the larger the coupling. 

2) At very small values of the coupling (p ~ 10 or greater) the fall is not felt at maximal 
rapidity y — 50 chosen for the evolution. This is in full accord with the predictions of earlier 
studies |H El E], in which it was concluded that the behaviour should be ~ exp(— AEy) 
with AE given by (jTDHjl . With p, = 1 and p = 10 one finds AE ~ 10~ 21 , so that the fall of 
the propagator should be felt only at fantastically large rapidities. 

3) However already with p, = 1 and A = 1/3 the decrease of the propagator with the 
growth of y is visible (see the non-logarithmic plot in Fig. 3). 

4) With small values of p, at A = 1 the propagator goes down practically as exp(— y), so 
that the triple pomeron coupling takes the role of the damping mechanism 

5) No singularity is visible at p, = and fixed A = 1. The curves with very small 
p, = ±10 -4 look completely identical. Of course this does not exclude discontinuities in 
higher derivatives, which we do not see numerically. 
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We next compared our exact results with two approximate estimates. The first is the 
propagator obtained in the approximation of 'big loops', Eq. ( IBB]) . The second is the 
asymptotic expression at small A obtained in earlier studies 

P(y) = p 2 e~y AE (119) 

with AE given by (11031) . The results for p = 1 and A = 1/10 and 1/3 are shown in Figs. 
2 and 3 respectively. As one observes, both approximations are not quite satisfactory. For 
quite small A = 1/10 the asymptotic limit, which we followed up to y ~ 50, is correctly 
reproduced in both approximations. However the exact propagator seems to reach this limit 
at still higher values of y (to start falling at fantastically large y ~ 10 20 ). The behaviour at 
small values of y is of course not described by f II 1 9 j) at all. The 'big loop' approximation 
works better and describes the rising part of the curve with an error of less than 10%, which 
goes down to 2% for y > 15. At larger values of A the situation becomes worse. As seen 
from Fig. 3 (A = 1/3) both approximations work very poorly. The 'big loop' approximation 
fails to reproduce both the magnitude of the propagator and its fall at large y. It describes 
the propagator more or less satisfactorily only at small values of y < 2. The approximation 
(11191) describes the trend of the propagator at higher y better but also fails to reproduce 
its magnitude. For still higher values of A both approximations do not work at all. So our 
conclusion is that both approximations are applicable only at quite small values of X/ p < 1/5 
and that the approximation of 'big loop' is better, since it also reproduces the growth of the 
propagator at smaller y from its initial value -P(O) = 1 to its asymptotic value p 2 at maximal 
rapidities considered. 

In Fig. 4 we compare the exact pomeron propagator with the calculations based on the 
perturbative expansion (up to two loops) and with the free propagator ( 11181) . One observes 
that there exists a region of intermediate values of py, from 1 to approximately 3 where, 
on the one hand, the influence of loops is already noticeable and, on the other hand, the 
perturbative approach works with a reasonable precision. 

7.2 hA amplitude 

In Fig. 5 we show our results for the amplitude A(y, g) at comparatively large rapidity y — 10 
as a function of coupling g to the target for different sets of p and A. This dependence mimics 
that of the realistic hA amplitude on the atomic number of the nucleus. In all cases except for 
(p, A) = (—1, 1/10) the amplitude A(g) rapidly grows as g rises from zero to approximately 
unity and then continues to grow very slowly at higher values of g clearly showing signs of 
saturation. To make this behaviour visble we use the log- log plot. In the exceptional case 
(p, A) = (—1, 1/10) one expects a similar saturation but at higher values of g, as will become 
clear in the following. 

We again compare our exact results with predictions of pure fan approximation, Eq. (1591) 
with gi = 1 and u = g, and the asymptotic formula similar to Eq. fll 19j) 

A(y,g)=p(l-e-^)e-y AE . (120) 

The results are shown in Figs. 6,7 and 8 for (p, A) = (1, 1/10), (1, 1/3) and (-1,1/10) respec- 
tively. With p positive the effect of loops is considerable, so that pure fan model gives a bad 
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description, especially at low values of g and not very small values of A. The asymptotic 
formula (120) works much better. At /x = 1 and A = 1/10 it gives very good results, as 
seen from Fig. 6. However at larger values of A//i its precision rapidly goes down. The 
case /i = — 1, A = 1/10 illustrated in Fig. 8 was chosen to see the influence of loops for 
the subcritical pomeron, where such influence should be minimal. In this case the pure fan 
formula gives the amplitude 

A(y = 10, g) = - r (121) 

vy i)> 1 + 100(1 - e~y) y=w K ' 

which grows linearly with g up to values around g = 10 and only then saturates at the value 
0.1e~ 10 . This explains a somewhat exceptional behaviour of A(y = 10, g) for such /i and A. 
As follows from Fig. 8 the influence of loops is in fact small. However it is greater than 
~ A 2 =l% as one could expect on simple estimates of a single loop contribution. In fact the 
pure fan formula overshoots the exact result by about 20%. 



8 Conclusions 

As mentioned, the zero-dimensional RFT was studied in detail some thirty years ago. Our 
aim was to clarify some mathematical aspects of the model and also to study it at different 
values of the triple pomeron coupling and not only at small ones. 

We have found that the model can be consistently formulated in terms of the quantum 
Hamiltonian. This formulation is valid both for negative and positive values of /i, in contrast 
to the functional approach which does not admit positive fi. The Hamiltonial formulation 
does not need to introduce a scalar product in the representation in which the creation 
operators u are diagonal. The wave functions seem to be not integrable in the whole complex 
u plane. Rather the standard scalar product should be used in which the creation and 
annihilation operators are expressed by Hermithean operators. 

A clear physical content of the theory is well seen after transformation to an Hermithean 
Hamiltonian made in [I]. However such transformation can only be done separately for 
positive and negative u, leading to two branches of the spectrum for the initial Hamiltonian, 
which include both positive anfd negative eigenvalues. The choice between two branches is 
determined by the signs of A and external coupling constants. 

For fi > a simple case of pure fan diagrams admits an explicit solution and reinter- 
pretation in terms of fashionable reaction-diffusion approach. Using this one can obtain an 
approximate expression for the pomeron propagator with loops of the largest extension in 
rapidity taken ito account, which gives a reasonable approximation at small values of A. 

Finally we performed numerical calculations of the amplitudes. They show that in all 
cases the triple pomeron interaction makes amplitudes fall at high rapidities. This fall starts 
later for smaller A and at very small A begins at asymptotically high rapidities (for A//i < 1/4 
it is noticeable only at fiy > 100). At small A the behaviour in y is well described by formulas 
obtained in earlier studies. However when A is not so small all approximate formulas work 
poorly. Our numerical calculations have also shown that there is no visible singularity at 
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[i = 0, in spite of the fact that the functional formulation meets with a divergence at this 
point. 

We have proposed a new method for the perturbative treatment of the theory, which may 
have appliciations to more sophisticated models in the world of more dimensions 
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Figure 1: The full pomeron propagator as a function of rapidity for different sets of \i and A 
indicated in the figure 
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Figure 2: The full pomeron propagator as a function of rapidity for ji = 1 and A = 0.1 
(the lower curve) as compared to the 'big loop' formula Eq. (88) (the middle curve) and 
asymptotic expression (119) (the upper curve) 
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Figure 3: The full pomeron propagator as a function of rapidity for /j, — 1 and A = 1/3 (the 
lower curve) as compared to the 'big loop' formula Eq. (88) (the upper curve on the right) 
and asymptotic expression (119) (the middle curve on the right) 
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Figure 4: The Pomeron Green's function as a function of \xy for /j, — 1 and A = 1/10: from 
bottom to top: the exact numerical evolution, the two loop perturbative approximation, the 
free case (e yfM ). 
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Figure 5: The hA amplitude at y = 10 as a function of the coupling to the tar- 
get. Curves from top to bottom on the right correspond to (/i, A) = (l, 0.1), (1,1/3), (1,1), (- 
1,0.1), (0.1,1), (0,1) 



25 



12 



10 - 




1 ' ' ' ' ' ' ' ' ' 1 

0.01 0.1 1 

g 

Figure 6: the hA amplitude A(y,g) at y — 10 for // = 1, A = 0.1 (the lower curve) 
compared to the pure fan prediction (the upper curve) and the asymptotic expression 
(120) (middle curve) 
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Figure 8: The hA amplitude A(y,g)ai y = 10 for /j, = — 1, A = 0.1 (the lower curve) as 
compared to the pure fan prediction (the upper curve) 
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